Stochastic limit approximation for rapidly decaying systems 
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The stochastic limit approximation method for "rapid" decay is presented, where the damping 
rate 7 is comparable to the system frequency Q, i.e., 7 ~ i7, whereas the usual stochastic limit 
approximation is applied only to the weak damping situation 7 ^ i7. The key formulas for rapid 
decay are very similar to those for weak damping, but the dynamics are quite different. From a 
microscopic Hamiltonian, the spin-boson model, a Bloch equation containing two independent time 
scales is derived. This is a useful method to extract the minimal dissipative dynamics at high 
temperature ksT S> TiO and the master equations obtained are of the Lindblad form unlike that of 
Caldeira and Leggett. The validity of the method is confirmed by comparing the master equation 
derived through this method with the exact one. 



One of the most attractive issues in fundamental 
physics is to understand dissipative dynamics from a 
microscopic point of view. This is not a trivial ques- 
tion since quantum theory is designed for closed sys- 
tems and has time-reversal symmetry. A standard ap- 
proach to this problem is to deal with the (dissipative) 
system of interest together with its surroundings, the "en- 
vironment," which has infinite degrees of freedom [Q-Q. 
The total system — "system" -|- "environment" — is then 
a closed system and is treated in a quantum-mechanical 
way. Caldeira and Leggett succeeded in formulating the 
"quantum Brownian motion" by making use of the path- 
integral method |^-^, and recent developments in tech- 
nology have revealed the validity of this approach in the 
field of quantum optics [|| . 

One of the powerful methods of deriving dissipative 
dynamics from the total system approach is the stochas- 
tic limit approximation method established by Accardi 
et al. |9|,[l0|. In quantum optics, for example, one is often 
interested in the decay of an excited atom in the radiation 
field. The total Hamiltonian reads 



Htot =Hs + \V + He 



(1) 



Hs and Hb being the Hamiltonians of the atom and the 
radiation field, respectively. Since the electromagnetic 
interaction W is weak (A <C 1), the weak-coupling limit 
A ^ gives a good approximation. Simultaneously mak- 
ing the coarse graining in time t — t j}? (A — > with r 
fixed), the stochastic limit approximation extracts the 
minimal (but notable) dissipative dynamics |pl]-p3[. 

However, the usual stochastic limit approximation ap- 
plied to the Hamiltonian (|^) gives a rather restricted dy- 
namics, i.e., weak damping. The lowest-order contribu- 
tion of the interaction to the damping rate 7 is propor- 
tional to A^ , and it is much smaller than the characteristic 



frequency of the system J7, i.e., 7 ^ f2, which is the case 
of quantum optics. For faster decay 7 ~ 12, which we 
call rapid decay, a separate treatment is needed. 

The rapidly decaying dynamics is also of physical in- 
terest, because it includes the case of quantum Brown- 
ian motion . It contains a richer variety of transient 
dynamics in contrast to the situation of weak damp- 
ing, where frequent oscillations of the system allow the 
"rotating-wave approximation," which simplifies the dy- 
namics considerably p^ ]. 

In this paper, we explore the stochastic limit approxi- 
mation for rapid decay. The total Hamiltonian is P,p^ 

Htot^X'Hs + XV + HB. (2) 

The system frequency f2 is, in turn, of order O(A^), and 
the resultant dissipative dynamics exhibits rapid decay 
7 ~ 12. It is shown that the key formulas are different 
from but similar to those obtained in the usual stochas- 
tic limit approximation for weak damping. This brings 
us with a Bloch equation with two independent time 
scales from a microscopic Hamiltonian, i.e., the spin- 
boson model. Furthermore, the master equations derived 
through this method are shown to be of the Lindblad 
form which ensures the conservation and positiv- 
ity of probability, unlike [p[ p^jl7| ] that of Caldeira and 
Leggett @. 

For the sake of simplicity, we shall focus on the fol- 
lowing spin-boson model [ p^ . The extension of the spin 
system to a general iV-level system is straightforward. 
The total Hamiltonian is given by Eq. (|) with 



Hp 



dkhujka\ak, (3a) 



V = ihaz / dk ( gkUk 



* t 



(3b) 
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The spin system Hs, which has two energy eigenstates ^^/(^) = J2 (^"^"(^) - -DlfolCr)) W/(r). (10) 



Hs\±) = ±\hQo\±), hf2o = Vs^ + A^, (4) 

interacts with the boson system Hb through the hncar 
couphng interaction V. The particular choice of the cou- 
pHng here is not essential except the linearity; gk is the 
coupling constant for the boson of mode fc, and through- 
out this paper it is assumed, as usual, that, before time 
t = 0, the two systems are uncorrelated and the boson 
system is in the thermal equilibrium state at tempera- 
ture T. The initial state of the total system /9tot(0) is 
thus given by 



Ptot(O) ^ PS® PB, 



(5a) 



PB 



where fcs is the Boltzmann constant. This model was 
discussed by the usual stochastic limit approximation in 
Ref. [0 but with the total Hamiltonian (Q). 

With Hq = X^Hs + Hb being the free part of the 
Hamiltonian, the time-evolution operator in the interac- 
tion picture, Ui(t), satisfies the Schrodinger equation 



d 
dt 



Ui{t) 



-\Vi{t)Ui{t), Ui{0) = 1, 



Vi{t) — giHot/hy^-iHot/h 

= zhJ2i^DMt)-DiAi{t) 



(6a) 



(6b) 



Here Da {a ~ {+, — ,0}) are the basic operators of the 
spin system, 

D± = \±){T\, Do = \+){+\-\-){-l (7) 
and Aa{t) the operators of the Bose field, 

Aa{t) = rffcgfeflfee-^^'^^-^''^")*, (8) 

with C+ - (+1-7,1-) =. A/hf2o = CI, Co = - 
— {—\o'z\—) = £/^^o, = ±^^0, and uq = 0. Coarse 
graining in time is made in the Schrodinger equation (|^) : 
Changing the time scale from the microscopic time t to 
the macroscopic one r = X^t, 



= Y. [d.\a,{tI\^) - D* i/lt (r/A=)) U,{rl>?). 

(9) 



Note that at the present stage, Eq. (|T^ should not be re- 
garded as mathematically fully justified. A more rigorous 
analysis requires a prescription to handle normal order- 
ing We shall see later that a consistent procedure 
of normal ordering can be obtained for the Heisenberg 
equations of the dressed spin system operators. 

Let us compute the correlation functions of the oper- 
ators of the Bose field 6q(t) in the thermal state pB- 
These operators play the role of the "quantum noises" in 
the macroscopic time scale. First notice the basic formula 

hm r drF{T) r |^G(a;)i,e-(-"^^^)-A^ 



lim / dtF(XH)G(t)e 
F(0) lim GiX^Cj). 



iX ujt 



(11) 



G{t) is the Fourier transform of G{lu). The correlation 
functions in the stochastic limit read 

{ba{r)bl{T'))B 

= hm r+„,(A2^„„,)e*^"°""°'^"^(^ - r'), (12a) 



{bl{r)ba'{r'))B 

= hm (r-„,(A^a;„„0)*e-^(— =')-<5(r - r') (12b) 

with the spectral functions 

rL'i^) = (l + ^(^)) ^a«'(u;), (13a) 



(13b) 



raa'{i^) = CaC:,r{Lu) = 2nCaCl, / dk \gk?5{u;k - w), 

(13c) 



the Bose-Einstein distribution function 

1 



N(uj) 



(14) 



we take the stochastic limit A — > and obtain 



and the shorthand notation lOaa' — (wq -I- L0a')/'2.. The 
correlation time of the operators of the Bose field 6q(t) 
is negligibly small in the macroscopic time, which makes 
it reasonable to call them quantum noises. 

The factor lim;^_»o r^a' {^^'^aa') is sensitive to the form 
of the spectral function Faa' (w) for small positive oj. [For 
w < 0, r'aa'(w) = by the definition in Eq. (13c).] Not- 
ing that A^(a;) ~ kBT/Tiu) for u ~ 0, one sees that the 
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only spectral functions that are linear (Ohmic) for small 
positive uj, 



(15) 



yield nontrivial values for \im\^o r^^,{X'^LOaa')- For 
ujaa' > 0, the zeroth order of the expansion in 

r^^>{X^uJa.a')/{kBT/n) = rj^a' + 0{X^nuJc.c.'/kBT), 



(16) 



survives in the stochastic limit A 



n 

2h 







0, so that 

{uJaa' > 0) 

{uJaa' = 0) 
(uJaa' < 0). 



(17) 



The case uJaa' = is treated separately since F^^, {to) is 
discontinuous at w = 0, and r^^'i^) defined here by 
^i'(O) = [ri'(0+) + /^i.(0-)]/2 for convention. 

The correlation functions ( p^ ) are similar to their coun- 
terparts in the usual treatment for weak damping |p|,p^, 
which are reproduced by the replacement X^uja 
in Eqs. (p^). There is, however, a significant difference 
between them: In the usual stochastic limit for weak 
damping, the rapidly oscillating factor e*('^°^'^a')'^/'^ ex- 
tracts the uja = '-^a' contribution only and makes the 
"rotating- wave approximation" exact |lj], while in the 
present case, the rotating- wave approximation is not ap- 
plicable in general. 

Now let us derive the Heisenberg equations for the 
dressed spin system operators in the bath, 



(18) 



from the Schrodinger equation ([lOl). The above set 
of equations yields the Bloch equation and is immedi- 
ately translated into the master equation. Differentiat- 
ing both sides of Eq. (^8|) with respect to r, and apply- 
ing the Schrodinger equation ([lo| ) to the time derivative 
of Ui{t), we will obtain the equations for 2?q(t) after 
computing the partial trace with the thermal state in 
Eq. (|l8|). Since we are working in the operator formalism, 
we rely upon the techniques of Thermo Field Dynamics 
(TFD) jl^ for the evaluation of the trace with the ther- 
mal state. In the language of TFD, the thermal state is 
represented by "thermal vacuum" \9). It is annihilated 
by the annihilation operator and its tilde conjugate 
^k, which are related to the operator by 



(19) 



Thus the key formulas for the derivation of the Heisen- 
berg equations for X'a(T) are the commutation relations 



between Ui{t) and the operators Xa{T), Xa(T), defined 
as the stochastic limit of the operators Sa{T/X?)/X and 
S^{t/X^)/X: 



Sa{t) = / dfc/TTiVK)5fcefce-^("'=-^'"=)*, (20a) 



i^{t) = C: / dfcv/iVK)g*|fee^("'=-^'"°)*. (20b) 



Notice that they are related to bair) by 5c((t) = Xq(t) -|- 
xUr). 

The commutation relations enable us to make normal 
ordering and to evaluate the partial trace. They are cal- 
culated in a way similar to that in the usual stochastic 
limit for weak damping |^,|l^: Observe the commuta- 
tion relations between S'q(t/A^)/A, Sa{T / X^) / X and the 
iterative solution 

Ui{t/X^) = 1 dr'jVi{r'/X')Ujir'/X% (21) 

and then take the stochastic limit. The key commutation 
relations are obtained as 



[Xo(r),iY,(r)] = - J] lim(^i:+„,(A2a;„„. 

a' 

X e*('^°-"°')^L'^,^Y7(T), (22a) 
[Xa(T),W/(r)] = ^ hm(zr-„,(A2a,„„0)* 

a' 

X e-*("°-""')"i?„.W/(T). (22b) 
Here S^^, (to) is the self-energy function, given by 

/•OO /'OO 1 [ 

= dt |^ri,(^')e-^(-'-)* (23a) 



io 27r ""^ 'u-uj'- 



^0^ 



(23b) 
(23c) 



whose imaginary part r^^, (w) will give the decay rate, 
and real part 



""^^ Jo 27r ""^ 'uj~u' 



(23d) 



will contribute to the energy shift in the following Bloch 
equation or in the master equation. The counterparts 
in the usual stochastic limit for weak damping |^,^ are 
again reproduced from Eqs. ( p2| ) by the same replacement 
X^uja — > as that for the correlation functions (E3) . 

By making use of the commutation relations (|22|), the 
Bloch equation for rapidly decaying spin system is ob- 
tained: 
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-if2a £^772 


^2^72 


\ 










-A^-fO 


eA-y'^ 




2^0 (r) 




V 




£^772 









(24) 



where e = e/hfio, A ~ A/fif2o, and 
7*^ = 2rjkBT/h. 



(25) 



Through the relations {±\ps{t)\^) = trs[ps2?zp(T)] and 
(±|ps(r)|±) = {1 ±trs[ps2?o(r)]}/2, this gives the mas- 
ter equation for the density operator of the spin system, 
Ps{t) = tr_B ptot(T), 

^psir) - -^[Hs,psir)] ^K, [a.,psir)]]. (26) 

Note that, for the model considered here, the energy shift 
disappears in the stochastic limit. 

The key formulas for weak damping and those for rapid 
decay are quite similar: The formulas for weak damping 
are reproduced by the formal replacement A^Wq. — > Wq. in 



those for rapid decay in Eqs. (12) and (p2|). However, in 
the case of weak damping where the total Hamiltonian 
is given by Eqs. (^ and (^), the Bloch equation is much 
simpler than that for rapid decay (|2^). The matrix in 
Eq. ( p^ ) for weak damping is diagonal (or decoupled), 
and furthermore, the damping coefhcients of 'D±{t) and 
VoiT) are given by — A'^Y/'^ and 7^ = A'^Y^ respec- 
tively, but with Y = r{f]n) coth{hf]n/2kBT) instead of 
that for rapid decay ( p5| ) [Q, i.e., the relation td = 2Tpi 
always holds between the decoherence time to = (7^ ) ~ ^ 
and the thermal relaxation time r/j = (7^)"^, while 
it does not always hold in spin-relaxation experiments. 
The phenomenological Bloch equation |2^ is given two 
independent time scales tji and to. The simplicity of 
the weak damping formulas is due to the rotating-wave 
approximation, arising from the frequent oscillations of 
the spin system 6**^'^°""°'^'^/''' , which suppresses many 
terms in the equations The Bloch equation (|2^ ) 

for rapidly decaying systems, on the other hand, con- 
tains a richer variety. Between the two decay constants, 
which are the real parts of the eigenvalues of the ma- 
trix in Eq. (p4|), there is no trivial relation as that for 
weak damping. Even more, there is a case, depending 
on the parameters, where three decay constants (three 
real eigenvalues) exist. Details of these features will be 
reported in elsewhere. 

It is easy to show that the master equation ( p6| ) has 
the unique thermal equilibrium state poq, which is pro- 
portional to 1: ps{t) — > pcq oc 1 as T — i- cx). pcq is 
nothing but the thermal state at infinitely high temper- 
ature. This is because the master equation (^6|) is valid 
for the situation X^hfio/kBT ^ 1. Remember the ex- 
pansion (p^. The stochastic limit approximation is a 



method to extract the minimal dissipative dynamics from 
the full of it. It picks only the processes where one boson 
is emitted or absorbed, and neglects the higher-order con- 
tributions . As for the stochastic limit approximation 
for rapid decay, one can say, in addition, that it extracts 
the dissipative dynamics at a much higher temperature 
than the characteristic energy scales of the system. 

The master equation ( p6| ) derived through the stochas- 
tic limit approximation for rapid decay is worthy of note. 
It is of the Lindblad form, which ensures the conserva- 
tion and positivity of probability [|5|. For the case of 
weak damping, many authors derived master equations 
of the Lindblad form, while for the case of rapid decay, 
much care is required The master equation de- 

rived by Caldeira and Leggett for example, is not of 
the Lindblad form |||,|l^,|rj . Such master equations lead, 
in certain cases, to unphysical results [pT| . 

It is interesting to look at the master equation for the 
model whose total Hamiltonian is given by Eq. (|) with 



1 

2M 



V 



-ih\j J dk {^guak ~ , 



(27a) 



(27b) 



and the same Hb as in Eqs. (pa|). The master equation 
derived through the stochastic limit approximation for 
rapid decay illustrated here reads 



rfr 



Ps{t) = --[hIps{t) 



MHorjkBT 
2h^ 



X, [x,ps{t)]], 
(28a) 



where Hg is the renormalized Hamiltonian with the 
renormalized frequency 



22 f° 
fi — .r?Q — 2f?Q I 

Jo 



duj r{uj) 
2t: lo 



(28b) 



The master equation ( 28a ) is of the Lindblad form since 
it lacks the term —{in{)rj/AK)[x, {p^ Ps{'^)Y\ which is con- 
tained in the one derived by Caldeira and Leggett 
This term may be neglected, when compared to the last 
term of the master equation (28a) if the temperature T 
is high enough. The same situation is found in the spin- 
boson model. The master equation derived by Munro 
and Gardiner in Ref . [ p^ , which is again not of the Lind- 
blad form, is reduced to the Lindblad form master equa- 
tion (pfl) in the high-temperature limit. 
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The model (|27|) treated above is exactly solv- 
able P,|6|,[7|Jl4|] , and it is possible to write down an exact 
master equation for the reduced density operator [Q . Let 
us finally confirm the validity of the master equation ( |2^ ) 
derived through the method presented here by compar- 
ing it with the exact one. The master equation, which is 
exact for an arbitrary A, is obtained through the method 
sketched in the Appendix and reads 



the reduced density matrix is given by (a;|ps(r)|a;') = 
J^oo'^^o J^^dx'oJ{x,x',T\xo,x'o,0){xo\ps\xQ)]^ is car- 
ried out without any approximation to give 



1 



J{Xf,x'j:,T\x^,x'i,0) = 



X exp 



J2 (^^kAki{T)Si - ^Ak0ki{r)Ai 



k,l 



, (Al) 



^f^^Psir) = [H',{T),psiT)]~tD,,{T)[x,[x,psiT)]] 

~ 2iDa;p{T)[x, [p,Ps{t)]] + r^p{T)[x, {p, ps{t)}], 

(29a) 



where (xk + x'^)/2, Ak = Xk - x'^. {k = {/,«}), 



yl//(r)=.l,,(T)=Af-ln|i^(r)|, 



(A2a) 



where 



(29b) 



The time-dependent coefficients f2ii{T), Dxx{t), Dxpir), 
and Fxpir) are given by Eqs. ([A^). They are reduced in 
the stochastic limit A ^ 0, as shown in the Appendix, to 



Dxpir) 0, Fxpir) 



0, 



(30) 



and the master equation ( p8[ ) is thus reproduced from 
the exact one (p9[), which shows the consistency of our 
method. 

We have discussed the stochastic limit approxima- 
tion for rapid decay and extended its applicability. Al- 
though some mathematical issues, e.g., the question of 
the convergence of the operators, remain unaddressed, 
the framework presented here is useful for practical cal- 
culations. 

The authors acknowledge useful and helpful dis- 
cussions with Professor I. Ohba. They also thank 
H. Nakazato and S. Pascazio for critical reading of the 
manuscript and enlightening comments, P. Facchi for 
fruitful discussions, and M. Miyamoto for advice. This 
work was supported by a Grant-in- Aid for JSPS Research 
Fellows. 



yl/.(r) = X^MF{T)^\n\F{T)l %(t) 



M 



X^F{t) ' 
(A2b) 



and 



r 



Om{t) ^ - I dri / dT2Ck{T,Tl)v{Tl ~ T2)ci{t,T2), 
^ Jo JQ 



(A3a) 



c;(r,r') ^ A^F(r') - ( ^ln|f^(r)| ) A^F(r'), (A3b) 



F(r)' 



(A3c) 



v{t)^MQo — F(X^Cj) coth^---^ cos Cot. (A3d) 
Jo 27r 2kBT 

F{t) is the solution of the integrodifferential equation 
F{t) + nlF{T) + 4t r dr'pir ~ t')F{t') = 0, (A4a) 



M 



p{t) = -MQq 



duj 
2^ 



F{X^Cj)sinCjT (A4b) 



APPENDIX: 

Here we briefly illustrate a derivation of the exact 
master equation ( |29| ) for the model (|^) with Hs^ V in 
Eqs. (E^), and Hb in Eqs. (pal), and show that the 



stochastic limit of it reproduces the master equation (28). 

The derivation of the exact master equation ( |29| ) 
is found in Ref. which is based on the influence 
functional approach proposed by Feynman and Ver- 
non The path integration of the time-evolution 
kernel for the reduced density matrix in the posi- 
tion representation, J{xf,x'j,Tf\xi,x[,Ti) [with which 



with the initial condition F(0) = 0, F{Q) — 1, and is 
given by an inverse Laplace transformation 



X^F(t) 



ds 



c 2m P + + 2fi(X^S)/M 



Jo 



(A5a) 



(A5b) 



C is the so-called Bromwich path. Note that time t has 
already been rescaled to r = X^t (and accordingly, the 
integration variable uj to uj — uj /X"^), but A is so far finite, 
and the solution (Al) is exact for an arbitrary A. 
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The exact master equation ( p9| ) is then derived by dif- 
ferentiating the kernel (Al) with respect to time r, re- 
placing Ai and Si with linear combinations of d/dAf, 
d/dSf, Af, and Ef multiplying by the initial density 
{xi\ps\x'^), and integrating over the initial coordinates Xi, 
x[. The coefficients of the master equation (|9|) are given 
by 



' M dT 



^..(r) = 4r,p(r) Off 



A. 



Aff 



(A6a) 



'7 



1^4,, (A6b) 



0^ 



and 



M 



(A6c) 



(A6d) 



Now we observe that the master equation (^8|) is re- 
produced from the exact one ( p9| ) in the stochastic limit 
A ^ 0. First notice that the function X'^F{t) is evaluated 
in the limit as 

A2F(r)-> / = ^ sin f2RT (A7) 



27ri §2 -I- nl 



with the renormalized frequency f^u given in Eq. (28b) 
Aki (t) are hence given by 



Ai 



Affir 

Af.ir), A, fir) ^ - 



MHr cot Urt, 
MQr 



(A8a) 
(ASb) 



and 6>fe((r) by 



sin i7i?T ' 

MQQrjkBT Qrt — sin i^^r cos (Irt 



h 



AQr sin Qrt 



MfiQ-qkBT Qrt cos Qrt — sin Qrt 



h 



inRsin^nRT 



(A9a) 



(A9b) 



It is easy to see that the substitution ( |A8D and ( [A9| ) 
brings us with those constants given in Eqs. and 
the exact master equation (|2^) is reduced to the master 
equation (psl). 
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